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Abstract 

We discuss the weak-coupling BCS theory of the heavy fermion supercon- 
ductor UPt3, accounting for that system's anisotropic, multisheeted Fermi 
surface by expanding the order parameter and pair potential in terms of ap- 
propriate basis functions of the irreducible representations of the crystal 
point group. Within a phenomenological model for the electronic structure 
of UPt3 chosen to capture the qualitative features of local density functional 
calculations and de Haas- van Alphen measurements, we show how Fermi sur- 
face anisotropy can favor pairing in certain symmetry classes, and influence 
the phase diagram of unconventional superconductors. We also calculate the 
Ginzburg-Landau coefficients, focussing on those coefficients relevant to cur- 
rent theories of the UPt3 phase diagram. 
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I. INTRODUCTION 



Several different scenarios have been proposed |T]-g to explain the unusual phase di- 
agram of the heavy fermion superconductor UPt3 in an external magnetic field |6]J7| and 
under pressure. This phase diagram consists of several superconducting phases of pre- 
sumably different symmetry. Almost all theories of the phase diagram involve unconven- 
tional superconductivity, i.e. superconducting order parameters of lower symmetry than the 
conventional one which only breaks the U(l) gauge symmetry. These scenarios have almost 
exclusively been discussed on the basis of generalized Ginzburg-Landau (GL) theories of su- 
perconductivity, where the free energy functional is derived purely by symmetry arguments. 
||lOHl3| It is the nature of this approach that such theories require a considerable number 
of adjustable phenomenological parameters. In addition the choice of the order parameter, 
which can be classified according to the irreducible representations of the normal state sym- 
metry group of the system, is open. Each representation corresponds to a Cooper pairing 
channel with a particular pairing energy, i.e. transition temperature. 

At present, no satisfactory microscopic theory exists which could assist in picking out the 
relevant representation and fixing the phenomenological parameters of the corresponding 
GL theory. Of course, it is possible to work based on semimicroscopic formulation of a 
generalized weak-coupling BCS theory. This type of approach is often used in the study of 
properties related to the energy spectrum of the quasiparticle excitations for unconventional 
superconducting states with a specific topology of nodal structure of the quasiparticle gap, 
and considerable insight for various low-temperature properties has been gained from these 



treatments. |T3[ However, they are insufficient to answer the open questions we have in 
determining a GL theory. One difficulty arises from the fact that the specific form of the 
attractive pairing potential is unknown. Various additional corrections to the simplified 
weak-coupling treatment might be important and need to be properly taken into account. 
One aspect is the shape of the Fermi surface, (FS) which has multiple sheets, some of them 
very anisotropic. ||I4|-|T6|1 In conventional weak-coupling calculations the FS has usually a 
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spherical symmetry. It would be helpful to connect the semimicroscopic and the GL approach 
by including effects of the FS anisotropy into the weak-coupling approach. Such effects are 
an important part in determining the symmetry of the relevant Cooper pairing channel as 
well as the values of the parameters in the GL theory, and have so far not been considered 
in the literature. 

Let us now briefly review the situation in the theory of UPt3. The scenarios for the 
UPt3-phase diagram can be divided into three classes. (1) The relevant order parameter 
belongs to a two-dimensional (2D) representation of the crystal point group symmetry D^h 
(hexagonal) of UPt3, or £?2(s,«)- 0,@|L7|i9 m this scenario a symmetry lowering field 

is required to lift the two-fold degeneracy of the 2D order parameter in order to describe 
the experimentally observed double transition. Such a field is given, for example, by the 
staggered moment of the weak antiferromagnetic order which appears below T/v ~ 5K 
in UPt3. |19| To this class of theories belongs in spirit an idea due to Zhitomirskii and 
Luk'yanchuk, |J which states that the system is actually close to an isotropic one. Hexagonal 
anisotropy is then taken to split the representations of the rotation group, leaving nearly 
degenerate ID or 2D representations. (2) The double transition can also occur for a triplet 
order parameter described by the one- dimensional (ID) representations of D 6 h group, when 
the spin-orbit coupling interaction is negligible. In this model the degeneracy due to the 
spin direction is lifted by the staggered magnetic moment. p],|2"0f (3) It is also plausible, that 
two irreducible representations of might have their transition temperatures very close 
to each other (accidental near degeneracy) such that they would generate the feature of the 
double transtion. |||5| The staggered moment does not play any role in this type of scenario. 

As a test case for these scenarios usually the tetracritical point in the if-T-phase diagram 
of UPt3 at (H,T) ~ (500Gauss, OAK) is analyzed. This point is most clearly observed for 
fields parallel to the basal plane, and is rather insensitive to the relative orientation with 
respect to the basal plane crystalline axes. The scenarios (1) have, in general, more difficulty 
than the scenarios (3) to reproduce this detail because of the anisotropy introduced by the 
staggered moment. One argument to overcome this flaw is the assumption that the staggered 



moment always follows the orientation of the applied magnetic field. A further problem for 
the type 1) theories is that the observed tetracritical point for fields along the c-axis is not 
well reproduced. Sauls suggests that the problem could also be solved if certain parameters 
in the GL theory are very small. |I] Based on this idea he gave an argument in favor of the 
irreducible representation E 2u (f-wave) which is a triplet pairing state. Joynt [^TJ argues 
similarly that in an E\ g representation an apparent tetracritical point for H || c is obtained 
if the system is nearly particle-hole symmetric. The data for this configuration are not 
conclusive. || The scenarios (2) of the order parameter given by ID representations and (3) 
which neglects the coupling between the staggered moment and the superconducting order 
parameter predict a tetracritical point in all directions. 

Strong evidence in favor of theories (types 1) & 2)) which rely on the antiferromagnetic 
order as a symmetry breaking field is the observed simultaneous disappearance of both 
the high temperature, low field "A phase" and the staggered moment with the application 
of hydrostatic pressure. p],|9| The disappearance of the two transitions could be taken as 
evidence for the Zhitomirskii-Luk'yanchuk scenario, suggesting a nearness to "isotropy". 
However the simultaneous vanishing of the staggered magnetization seems more natural in 
the context of the 2D theories based on a symmetry breaking field. More recently, it was 
found that at the pressure value where the transitions merge, application of basal plane stress 
"resplit" the transition, lending further support to these scenarios. |22|] On the other hand, 
in these theories the staggered magnetization generally must be assumed to rotate with the 
applied field in order to understand both the isotropy of the phase diagram for H _L c and 



the weak sixfold variation of the upper critical field H c2 with angle in the plane. |S| The 
recent experiment of Lussier et al. which found no rotation of the staggered moment with 
field appears to contradict this hypothesis, however. [p^| 

Further constraints on the theory are provided by recent Knight shift measurements, 
which show no temperature dependence for any direction of the field below T c . [^] This 
striking result would be consistent with an equal spin pairing state such as in 3 He — A, with 
Cooper pair spins free to rotate with the external field. Such a model (type 3), in which 



the order parameter is a spin triplet state with 3 spin components split by the staggered 
moment, has indeed been proposed by Machida and co-workers. 0,^6j This theory has the 
disadvantage of generically exhibiting three transitions in zero field rather than two, due 
to the three independent order parameter components. It also cannot explain the apparent 
Pauli limiting of the upper critical field for H || c. [E7 



Thus, the existing experimental data have not been particularly kind to any of the 
current theoretical scenarios. It is not our aim to identify the correct theory here. Rather, 
we focus on an aspect which has been generally neglected in oversimplified models, namely 
which types of order parameter might be favored by the unusual, highly anisotropic Fermi 
surface of UPt^. Within an extension of the semimicroscopic weak-coupling theory, we 
hope to be able to suggest which scenario is most likely, although no concrete proof can 
be given. From our analysis we expect to answer or shed light on some of the questions of 
interest. Are two nearly degenerate order parameters very likely to appear in UPt3? Does 
the scenario (1) given by Sauls satisfy the conditions that GL theory would produce the 
isotropic tetracritical point? Is the system close to "isotropy" in the sense of Zhitomirskii 
and Luk'yanchuk? Which representations are favored by hexagonal anisotropy? And, how is 
the pairing potential distributed among the various sheets of the multisheeted Fermi surface? 

Our discussion is considerably simplified by the assumption that a weak-coupling type 
of theory can be applied in case of UPt3. Several properties of UPt3 support the valid- 
ity of such an approach. First, UPt3 (in contrast to UBei3) is in a rather well-developed 
Fermi liquid (FL) regime at the onset of superconductivity. For example, the linear specific 
heat and the Independence of the resistivity are suggestive for the existence of well-defined 
(Landau) quasiparticles. Their effective mass is extremely large such that the de- 



generacy temperature Tp is of the order of 10 ~ 2§K. The question arises whether strong 
coupling effects are important in such a situation. They are measured by the ratio T c /Tp 
which is here of the order of 1/20 . Thus we neglect these corrections, although they are 
considerably stronger than in conventional superconductors with T c /Tp ~ 10~ 3 . Indirect 
justification for this assumption is provided by the size of the specific heat discontinuity at 



the upper transition point, AC/Cn ~ 0.8, which has to be compared with the BCS 
weak-coupling result, AC/Cn ~ 1.43. The FL renormalizations are of course quite large, as 
evidenced by the large effective mass. However, we will focus our attention to the proper- 
ties of the superconductor close to the transition temperature, where these corrections are 
irrelevant. [|3D| Furthermore we simplify our discussion by taking the pairing potential as the 
only interaction among the quasiparticles into account. Thus the only modification to the 
usual weak-coupling theories is inclusion of the Fermi surface anisotropy. This implies also 
that pair correlation can be confined to a narrow energy shell about the FS of a width u c , 
the cutoff energy associated with the electron-electron interaction (analogous to the Debye 
frequency u>d in conventional superconductors). 



II. FERMI SURFACE ANISOTROPY AND THE CHOICE OF THE ORDER 

PARAMETER 

A. The linearized gap equation 

The BCS description of the superconductor is based on the following typical Hamiltonian 

k,s Z k,k' si,..,s 4 

where c ks (c ks ) denotes the creation (annihilation) operator of a quasiparticle with momen- 
tum k and spin s. The quasiparticle band e k is measured relative to the Fermi energy Ep 
and V SlS2S:jS4 (k, k') is the pairing potential which is finite only for momenta k and k' for 
which |e k | < u c , < u c , where only scattering processes between particles of opposite 
momentum are included. The superconducting state is expressed by the mean field 

A ksiS2 = ^lS2S3S4(X k ')(Ck's3 C -k' S 4> ( 2 ) 

k' s 3 ,S4 

where (....) is the thermal expectation value. The mean field A kss / describes the gap of the 
quasiparticle spectrum E k = [e k + tr(A k A k )/2] 1 / 2 . The gap function A k is a 2 x 2-matrix 
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in spin space. By inserting the expectation value for a given temperature T in Eq. (2) we 
obtain the self-consistency equation for the gap function, 

A ksiS2 = -EE V SlS2S3SA (K k')%^tanh (_^L.) (3) 

This equation can be used to determine T c by considering the limit of the order parameter 
amplitude A — > which yields the linearized gap equation 

uA ksiS2 = v sis 2 s 3 s 4 (K k')A k/S[iS4 5(e k ) (4) 

where 1/v = ln(1.14o; c /A; B T c ) is obtained by integrating tanh(e / '2k bT c ) /2e over e in the 
energy shell, — u c < e < u c . This equation is an eigenvalue equation where the largest 
among the eigenvalues v determines the superconducting transition temperature T c and the 
corresponding symmetry of gap function. 

The Cooper pairing states can be classified into even (spin singlet) and odd (spin triplet) 
parity states which can be parametrized by a scalar function 

A k = A^(k)ia y (5) 

for even parity and a vector function 

A k = A d^(k)ia^a y (6) 

H=x,y,z 

for odd parity. If we assume that in the scattering described by V^ lS2SsS4 (k, k') spins are not 
affected then the pairing potential has the following form 

Vsis2S3S4 (k) k ) = V(k, k )(S Sl s 3 5 S2S4 i 5 SlS4 5 s , 2 s 3 ) (7) 

where the + (— )-sign occurs for V SlS2S3S4 (k, k') supporting odd (even) parity pairing. In this 
case the linearized gap equation can be expressed in the new parametrization as 

i^(k) = -£V(k,k>(k')(^kO (8) 
k' 

and 

<(k) = - £ V(k, k')d,(k')5(e k ,) (9) 
k' 

Note that in the case of odd parity pairing the different components of c? M (k) do not mix in 
the linearized gap equation. 
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B. Order parameter 



In our weak-coupling approach we ignored quasiparticles further away from the FS so 
that in the linearized gap equation given in Eq. (4) the momentum sum is restricted to the 
Fermi surface. In this case we express the gap function in terms of the linear combinations 
of the spherical harmonics defined on an appropriate Fermi surface. It is however necessary 
to choose this set of polynomials {^(k)} so that they are orthonormal to each other in the 
sense that the scalar product is determined as a FS average 

E^(k)^'(k)5(e k )=^ (10) 

k 

where / is a label for the polynomial. This property is most easily achieved by classifying 
them with respect to the irreducible representations of the crystal field point group, D 6 h 
in our case. Then only the polynomials within each irreducible representation have to be 
orthogonalized. 

The gap function for the singlet state can now be written as 

n 

^(k) = J>Mk) (11) 
i=i 

where rji are the complex expansion coefficients. For the triplet state we have to express the 
vector d(k) which requires an additional set of polynomials {</?f(k)} which leads to 

^(k)=Ewf(k) (12) 
i=i 

Orthogonality in this case is defined by 

E E ^(k)^(k)5( £k ) = 5 lv (13) 

k fj,=x,y,z 

In Tabs. I-III we catalogue these polynomials up to the order 3. We can now use them 
to express the linearized gap equation in matrix form as we show here for the case of even 
parity 

vTH = -'£{l\V\l')Th' (14) 
v 
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where 



(l\V\l>) = ^6fa)6(ev)<fi l (k)V(M,\c!)MV) (15) 
k,k' 

Clearly, this eigenvalue equation decays into block form where each block belongs to one 
irreducible representation of D^. Thus we may consider the problem for each representa- 
tion separately. There are representations in Tabs. I and III which have more than one 
polynomial. It is necessary to include all of them and ensure that they are orthonormalized. 
Note that the orthonormalization depends on the properties of the FS. 

Sometimes it is advantageous to express the gap function in terms of so-called FS har- 



monics, as introduced by Butler and Allen, |31| , |32|| and discussed in the context of the UPt^ 
pairing problem in Ref. ||33|| . The Fermi surface harmonics are homogeneous polynomials 
of the Fermi velocity, Vk = Vk^k and in the case of D 6 h point group they can be obtained 
from Tabs. I-III by replacing k (ki) by Vk (f j). These functions are to be orthonormalized 
in the same way as the polynomials constructed from spherical harmonics. We will use the 
FS harmonic representation of the order parameter in Section 3 of this paper. 



C. Pairing interaction and electron band 

Now we have to specify the pairing interaction V^k, k') and the electron band e^. In a 
spherical symmetric system the pairing states would be classified according to their angular 
momentum quantum numbers, I and m (and radial quantum numbers among which we 
will consider here only the lowest ones in each angular momentum channel). States for a 
given I are degenerate, i.e. yield in the linearized gap equation above the same transition 
temperature T c . It is natural then to write the corresponding pair potential as 

V SlS2S3S4 (k,k>) = E Y em (k)Y^(k')[6 slS3 5 S2Si - (-l)% 1Si S S2S3 ] (16) 

l m=-i 

where Y^ m are the spherical harmonics depending on the Fermi vectors k. Note that the 
Yi m (k) are dimensionless but are different from the usual Y em (k) away from the isotropic 
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Fermi surface point. In this notation Vi is negative for an attractive potential. In the 
following we will consider the FS anisotropy (crystal field) effects on the eigenvalues of Eq. 
(4). The information we can obtain is how the degenerate eigenstates evolve as the FS 
gradually becomes anisotropic. We will keep the interaction in this form where Vi is an 
unknown parameter. 

Information from experimental data as well as band structure calculation can help us 
to obtain an approximate electron band e^. While the FS of UPt3 is highly complex and 
consists of at least five sheets, []14|-|T6[1 not all of them are equally important. We will focus 
here on two bands only which are supposed to be the ones with the largest and the third- 
largest contribution to the density of states, that is the T 3 and T 2 FS sheet respectively 
(notation from Ref. 16). A simple but crude fit of the band structure data can be achieved 
by the following form for 

e k = a x k\ + a 2 (k 2 x + k 2 y ) + a 3 (k 3 y - 3k 2 x k y ) 2 - e F (17) 

where aj and Bp are fit parameters. The T2 band is fitted by an ellipsoidal which is rota- 
tionally symmetric about the c-axis, that is a 3 = in Eq. (|T7D. The hexagonal anisotropy 
is larger for the T 3 band which we approximate by a "wrinkled" ellipsoid with a 3 ^ 0. Note, 
that the latter band has additional pockets at Brioullin zone boundary which are not repro- 
duced by £k from Eq. (0). However, this form of has the advantage that we can follow 
continuously the development of the anisotropic bands starting from a spherical symmetric 
one (ai = a 2 ^ and a 3 = 0). The values of the parameters in for the anisotropic 



band are obtained using the ratios kp x /kp z and kp y /kp z from LDA calculations |16[ and by 
conserving the area of cross section in the k v -k z plane as observed in de Haas-van Alphen 



measurements. |T3| Our best fit parameters are given in Tab. IV, and the FS cross sections 
by the symmetry planes are shown in Fig. 1. This fit leaves still the density of states inte- 
grated over the FS as an undetermined parameter. In the following we will express various 
quantities in dimensionless units refering to the integrated density of states. 
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D. Solution of the linearized gap equation 



We address now the problem of the stability of the possible pairing states in the given 
T 2 - and T 3 -band. In reality, superconductivity is a combined effect of all FS sheets, because 
quasiparticle scattering among the different sheets will couple their pairing amplitudes. It 
would be not difficult to introduce such interband coupling in the pairing interaction V"(k, k'), 
but this would lead to additional unknown coupling parameters. In the following discussion 
we will avoid this complication and consider the effect of the two FS sheets separately as if 
they were uncoupled. 

We investigate the effect of the FS anisotropy on the pairing states by continuously dis- 
torting a spherical FS into the above defined FS sheets r 2 or (Eq. (|17|)). For the spherical 
FS all pairing states of the same angular momentum £ are degenerate. The degeneracy of 
these 2£ + 1 states is lifted when the FS deviates from the spherical shape, and they can 
be classified by the irreducible representations of the new, lower symmetry _D 6 /i- Restricting 
our attention to only one value of £ (£ — 1, 2 and 3) at a time we need only to take one 
coupling parameter, Vg, into account. This assumption may be too simplistic for realistic 
pair potentials, but we make it in the first analysis to isolate the effects of the Fermi surface 
anisotropy. We use an eigenvalue of Eq. (|4j) on the spherical FS of the same integrated 
density of states, N(0), as a reference value for v. Expressed in dimensionless units v' = vjvn 
does not depend on iV(0), which is not determined within this approach. The distortion 
from a spherical FS that is the effect of a crystal field is included by the dimensionless linear 
parameter t changing from to 1 , with t = corresponding to the spherical Fermi surface of 



the radius adjusted to dHvA data |14| and t = 1 representing T 3 or T 2 Fermi sheet (Tab. IV). 
Therefore an increase in t implies a higher hexagonal ^3) or ellipsoidal (r 2 ) perturbation of 
the Fermi surface. The t-dependence of the energy band coefficients for the hexagonal FS is 
given by 
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ai (*) / (e F c 2 ) = 0.225 - 0.068t 

a 2 (t)/(e F c 2 ) = 0.225 + 0.016t (18) 
a 3 (t) I (e F c 6 ) = 0.016t 

and for the ellipsoidal Fermi sheet 

a 2 (t) I {e F c 2 ) = 0.666 - 0.200t 

a 2 (t) I (e f c 2 ) = 0.666 + 0.285t (19) 
a 3 (t) = 

where c is the z-axis lattice constant (Tab. IV). The solutions of the linearized gap equation 
(Eq. (4)) for the p-wave, d-wave and f-wave pair potentials as the functions of the crystal field 
t are shown in Fig. 2 for the case of pairing interaction on T 3 Fermi sheet only (hexagonal 
crystal field) and Fig. 3 for the case of pairing potential on T 2 Fermi sheet (ellipsoidal crystal 
field). We present the critical temperature T c for triplet states (Figs. 2a, 2c, 3a, 3c) with d || z 
(Eq. (|(J)) only as the states with d±z are degenerate with them according to Tabs. I and III. 
Particularly interesting is the result for the d-wave states with the pairing interaction on T 3 
FS (Fig. 2b). We observe here a competition between A lg and E lg states as the hexagonal 
anisotropy of the Fermi sheet is increased. Finally, for our best fit T 3 FS (t = 1) the transition 
temperature of A\ g state (T c (Ai g )) is higher than that of E\ g state (T c (Ei g )), but the system 
is very close to the degeneracy point (t w 0.85) where T c (A lg ) = T c (E lg ). Thus application of 
a hydrodynamic pressure which decreases the crystal field t and restores "isotropy" may drive 
the system to a state with a single transition temperature T c = T c (Ai g ) = T c (E lg ), consistent 
with the Zhitomirskii and Luk'yanchuk scenario. The role of the hexagonal crystal field is 
worth mentioning, since as we show in Fig. 3b the cylindrical crystal field lifts the spherical 
FS (t = 0) degeneracy of d-wave states and splits the transition temperatures corresponding 
to different irreducible representations. Therefore the situation of near degenerate states 
Ai g and E\ g is not realized when the pairing mainly takes place on the ellipsoidal T 2 FS. 
Analysis of the f-wave states shows that the E 2u irreducible representation, suggested to 



describe the superconducting state in Sauls' scenario, |IJ,|Ij||34| is the one with the highest 
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transition temperature on T 3 FS but almost degenerate with A\ u representation (Fig. 2c). 
On the other hand, the transition temperature of A\ u state is significantly higher than that 
of E 2u on the T 2 sheet of FS (Fig. 3c). 

III. GINZBURG-LANDAU COEFFICIENTS 

Direct calculation of the GL coefficients incorporating aspects of the anisotropy of the FS 
is important for several reasons. The question of most current interest is the consistency of 
the Sauls E 2u scenario for the UPt3 phase diagram. |l| This theory relies on the smallness of 
a particular term in the GL gradient free energy, which may be shown to destroy the isotropy 
of the tetracritical point in the 2D scenario. In the particular case of the E 2u representation, 
this term may be shown to vanish in cylindrical symmetry, and it is therefore hexagonal 
anisotropy alone which is responsible for the disagreement between the predictions of the 
E 2u scenario for the phase diagram and experiment. If the effect of the anisotropy could 
be shown to be quite small, i.e. the system were in this sense quite close to a cylindrically 
symmetric one, one would have considerably more confidence in the theory. 

Before discussing this point in some detail, we remark however that other important 
physical information can be gleaned from a weak coupling calculation of the GL coefficients 
over the FS. Experimentally measurable quantities at the transition depend directly on com- 
binations of the GL coefficients, e.g. for the specific heat jump, AC/Cn = o> 2 /2f3\T c . If 
sufficient independent measurements can be made to pin down the coefficients individu- 
ally, and if the correct representation is known, comparison with weak-coupling calculations 
including the correct Fermi surface could in principle allow one to determine the size of 
strong-coupling corrections. This might in turn serve as a guide for theories of the UPt3 
normal state. 

The form of the GL free energy density has been given in several places, e.g. Refs. 1,13. 
For the 2D representations in hexagonal symmetry it takes the form 
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/ = a\rj\ 2 + /?i|r/| 4 + folwl 2 + 7iM 6 + l2\v\ 2 \ r l r l\ 2 + 

(20) 

Ei,j= x ,y{^ (PiVj) (diVj)* + k 2 (difji) (djrjj)* + k 3 {d^) (djrji)* + k 4 (d z rjj) (d z rjj)*} 
with the order parameter given by the 2D vector 

V = (Vi,V2) (21) 

that is for 

i) (k) = (k) + 7/202 (k) (22) 

in the case of even parity (Eq. (0)) and 

d, (k) = m tf (k) + rttf (k) (23) 

for an odd parity (Eq. @), where (f>i (k) and (pf (k) are the polynomials listed in Tabs. I- 1 1 1 . 
In the spin-triplet states we have assumed that the d vector is real and oriented along the 
crystal direction c (i-axis) by strong spin-orbit coupling in accordance with Sauls' scenario. 
|U This conjecture leads to the identical forms of the GL functionals for both even and odd 
states. 

Expressions for the GL coefficients are straightforward to obtain through expansion of 
the Luttinger-Ward functional for the BCS superconductor in powers of the order parameter, 
and identification of the prefactors of the various invariants present in the phenomenological 
form in Eq. (18). One finds for the case of even parity 
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T 

a = N In — 

J- r 



ft = 2(3 2 = 2(3 BCS (#(k)$00, 
n = ll2 = 2 lBCS (#(k)#(k)) 
«i = *»c(«?#(k)) (24) 
k 2 = k 3 = k wc (k) 4>j (k)) 

^4 = «wc (Vz4>i (k)) 
^123 = «W {V 2 i<t>1 (k)) 

where it is to be understood that i ^ j = x, y), K123 = «a + k 2 + k 3 , and the constants 
are given by 

pBCS = T6 as) J^f 

31 . . JV 

Ibcs = -^C (5) ^1 (25) 

For odd parity these expression are analogous where we have to take into account that a 
product is replaced by a scalar product (pi ■ <pj. 

Knowledge of the Fermi surface and of the basis functions for the various representations 
now allows one to evaluate, in principle, all coefficients of interest (Eq. fl2"4|)). There are 
several hidden problems in this analysis. The most serious has been alluded to above, 
namely our lack of knowledge of the relative pairing weight on each of the sheets of the 
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Fermi surface. As before, we approach this difficulty by assuming that pairing takes place 
either on the T 3 or the T 2 sheets of the Fermi surface, and discuss both cases separately. The 
second has to do with the choice of basis functions. We compare two different methods and 
calculate the coefficients for the expansion of the order parameter in the spherical harmonics 
and in the FS harmonics. Another complication arises from an infinite number of functions 
belonging to each irreducible representation in a crystal. Thus our choice of the <pi(k) and 
(/9/(k) is perforce somewhat arbitrary. In this work we generally work with the lowest order 
polynomial functions for each representation. We can test to see whether the choice of 
higher order polynomials will greatly alter the size of the coefficients in Eq. ( f24|) (it does 
not), but a more quantitative approach would allow the order parameter to adjust to 
the given angle-resolved density of states by taking a general linear combination of many 
such functions. This would have, for the FS harmonic representation particularly, the virtue 
of allowing for smoothing of singularities which inevitably result from the rapid changes of 
Vk over a complicated Fermi surface. Working with a single function (^(v^) and (pf (vjj is 
sufficient for our purposes, however, since the coefficients in Eq. Q23| ) are always integrated 
over the full FS under consideration, and only the degree of hexagonal anisotropy should be 
important. 

A special role in the topology of the H — T phase diagram of the the 2D representations 
is played by the coefficient k 2 (Eqs. (^0|), (|24]) ) which mixes the gradients of the two 
components of the order parameter. The solution of the eigenvalue problem determining 
the critical field lines in the H — T plane is analogous to a Schodinger equation for a charge 
moving in a magnetic field. The existence of a tetracritical point has been shown therefore 
to correspond to a level crossing, which cannot exist in the absence of a conserved quantum 
number. The k 2 mixing terms correspond to a level repulsion (or hybridisation) term in this 
analogy, and thus prevent a crossing of the critical field lines for H || c. 

By examining Table III, it is easy to see from Eq. (22) that the coefficient k 2 in a 
cylindrically symmetric system, where v || k, vanishes for E 2u symmetry whereas for the Ei 
representations, k 2 / k W cVf± is of Oil). Similarly it is easy to check that Ki/n wc Vp ± is of 0(1) 
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for both representations. Sauls' argument |1] is therefore that if hexagonal anisotropy effects 
are weak (k 2 <C the critical fields may come so close as to mimic an apparent tetracritical 
point, even for H || c. To test this proposition, it is interesting to know not simply whether 
K2 is in fact much smaller than Ki when evaluated over the UPt3 Fermi surface, but whether 
or not this result is "accidental" or due to the near cylindrical symmetry of the Fermi 
surface. Is the change in k 2 due to hexagonal anisotropy in fact small? Sauls proposes 
that it is, and gives a perturbative analysis for a band structure of type given by Eq. (|Ti 



assuming 03 small. In this case one finds k 2 /ki — 0.1 . |36| It is however not clear that 
a perturbative analysis is applicable. Here we calculate Ki, /t 2 numerically for the band 
given by Eq. (|17D, and display the results as a function of the anisotropy parameter a 3 in 
Fig. 4. In the calculation presented in Fig. 4a the order parameter was expressed in the 
spherical harmonics, whereas Fig. 4b shows the result for the FS harmonic representation 
of the gap function. The T3 Fermi sheet is given by a 3 /{eFC Q ) = 0.016 (Table. IV) and the 
corresponding Ki, k 2 values are marked with a dotted line in Fig. 4. We read from Fig. 4a 
that ^2/^1 — 0.46 for the order parameter expressed by the spherical harmonics and from 
Fig. 4b that 1^2/^1 — 0.73 for the FS harmonic expansion of the order parameter. Because of 
the cylindrical symmetry of T 2 FS, the k 2 coefficient calculated over this FS sheet vanishes. 



IV. CONCLUSIONS 

We have tried to take crude features of the UPt% Fermi surface into account within 
a weak-coupling theory of unconventional superconductivity to see which order parameter 
symmetries are favored by Fermi surface structure alone. In the current philosophy, these are 
classified according to their angular momentum quantum numbers in the fictitious spherical 
system obtained when the ellipsoidal and hexagonal deformations of the "true" Fermi surface 
are "turned off" . Within this scheme, we have examined states of p, and d, and / symmetry. 
Other important elements of the physics, such as the actual nature of the spin-orbit coupling 
and the range of the pairing, may also play major roles, but have been neglected here. The 
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lack of knowledge of a microscopic pairing mechanism-and consequent inability to specify 
the relative pairing weights on the various Fermi surface sheets-prevents us from drawing 
firm conclusions from an analysis of this type, but our results are suggestive. 

It is worth observing that the few representations we have found to be stabilized by 
Fermi surface anisotropy are precisely those under active consideration as candidates for 
the UPt 3 order parameter: the two-dimensional representations Ei g and E 2u and the ID 
representations Ai u and A lg . Both of the ID representations found are of interest in the 
accidentally degenerate representation scenario of Garg and co-workers, which requires, 
however, mixing of nearly degenerate v4-type and B-type representations. All of the latter 
are disfavored by the Fermi surface deformations we consider; our results do not therefore 
lend support to this model. The trivial representation A\ g corresponds to basis functions 
with the full symmetry of the Fermi surface, which typically are fully gapped. Such functions 
will be strongly supressed by local Coulomb interactions, which we do not consider here; 
they are by themselves ruled out by experiment. Exceptions occur if the system condenses 
in an A\ g state with nodes, or mixes with another representation supporting basis functions 
with nodes, as in the Zhitomirskii- Luk'yanchuk scenario. || The near degeneracy found in 
Figure 2b) for the A\ g and E\ g states for t ~ 0.85, which is very close to the best fit to the 
r 3 Fermi surface sheet (t — 1), is some support for this picture. 

We have attempted to account for our ignorance of the true distribution of the pair 
weights over the various Fermi surface sheets by examining models with pairing on the two 
sheets with highest density of states according to de Haas- van Alphen measurements, namely 
T 2 and T 3 . As T 2 is nearly ellipsoidal but T 3 has strong hexagonal deformations, this choice 
has the additional virtue of crudely separating hexagonal and ellipsoidal variations. This 
may be of significance for current theories: for example, the hexagonal deformation appears 
to be important in stabilizing the E 2u representation. Confirmation of the E 2u model would 
therefore be indirect evidence for pairing primarily on the T 3 sheet. Such conclusions could 
be used as guides for microscopic pairing theories. 
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The scenario based on the E 2u representation due to Sauls [JJ is based on the assumption 
that the level repulsion term in the Schrodinger-like equation for the critical field lines van- 
ishes due to the proximity of the system to cylindrical symmetry. This term is proportional 
to the ratio of GL gradient coefficients (^/fti) 2 , [[3(J which vanishes in this limit for the E 2u 
state. We have calculated these coefficients over our model Fermi surfaces fit to LDA cal- 
culations and dHvA measurements, and find that this term is not small, but varies between 
.21 and .53, depending on our assumptions regarding the exact form of the basis functions. 
While the parametrization of the Fermi surface we have adopted is extremely simplistic, it 
has the virtue that one can discuss the proximity to cylindrical symmetry. What Figure 4 
suggests is that the UPt 3 system may not really be regarded as close to cylindrically sym- 
metric, and that therefore if calculations of the GL coefficients over the true Fermi surface 



produce a smaller value, |36[ it should be regarded as coincidence. This undermines some- 
what the attractiveness of the Sauls E 2u scenario, but of course does not constitute proof 
against it. One way in which this conclusion could be avoided is to have contributions to k 2 
which change sign on different Fermi surface sheets, leading to a smaller effective k 2 . This 
effect is indeed found in a calculation with the full UPt 3 LDA Fermi surface and a model 



with weight distributed equally over all sheets. 
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TABLES 

TABLE I. Basis functions for the gap function: p-wave 
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TABLE II. Basis functions for the gap function: d-wave 
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TABLE III. Basis functions for the gap function: f-wave 
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TABLE IV. The functions and the coefficients used in fitting Y2 and T3 energy bands, 
c = 4.9027A is the z-axis lattice constant (Ref. 35). 
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a,i 


r 2 


r 3 
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ai/{e F c 2 ) 
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a 2 /(e F c 2 ) 
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FIGURE CAPTIONS 



Fig. I. T 3 (solid line) and T 2 (dashed line) Fermi sheet plot in the symmetry plane: a) 
ATM, b) LTM, c) HYK. 

Fig. 2. The hexagonal crystal field t effect on the critical temperature of the a) p-wave 
states, b) d-wave states, c) f-wave states, t — corresponds to isotropic, t — 1 to best fit 
dHvA-LDA T 3 Fermi surface. 

Fig. 3. The cylindrical crystal field t effect on the critical temperature of the a) p-wave 
states, b) d-wave states, c) f-wave states, t — corresponds to isotropic, t — 1 to best fit 
dHvA-LDA T 2 Fermi surface. 

Fig. 4. Normalized K\ (solid line) and k 2 (dashed line) stiffness coefficients as the func- 
tions of the normalized FS hexagonal anisotropy parameter 03: a) spherical harmonics rep- 
resentation, b) FS harmonics representation. The dotted line corresponds to the value 
G3A-FC 6 = 0.016 obtained by a fit to the T 3 Fermi surface sheet. 
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